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Abstract 

The magnetic relaxation processes in disordered two-dimensional ensembles of dipole-coupled 
magnetic nanoparticles are theoretically investigated by performing numerical simulations. The 
energy landscape of the system is explored by determining saddle points, adjacent local minima, 
energy barriers, and the associated minimum energy paths (MEPs) as functions of the structural 
disorder and particle density. The changes in the magnetic order of the nanostructure along the 
MEPs connecting adjacent minima are analyzed from a local perspective. In particular, we 
determine the extension of the correlated region where the directions of the particle magnetic 
moments vary significantly. It is shown that with increasing degree of disorder the magnetic 
correlation range decreases, i.e., the elementary relaxation processes become more localized. The 
distribution of the energy barriers, and their relation to the changes in the magnetic configura- 
tions are quantified. Finally, some implications for the long-time magnetic relaxation dynamics 
of nanostructures are discussed. 

PACS numbers: 75.75.+a, 75.60.Jk, 75.50.Lk, 
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I. INTRODUCTION 



Two-dimensional (2D) ensembles of interacting magnetic nanoparticles are currently 
the subject of an intense fundamental and technological research activity.-i^ The magnetic 
particles in these systems are subject to both local couplings and long-range interactions. 
By varying the interparticle distance the relative importance of these two contributions 
can be tuned to a certain extent. The magnetic response can then be correlated with 
the size distribution and the geometrical arrangement of the nanoparticles, the latter 
being experimentally accessible through a number of microscopy techniques. Of par- 
ticular interest is the slow magnetic dynamics,- like magnetic ageing, rejuvenation, or 
magnetic viscosity.- These non-equilibrium effects are reminiscent of the spin glass be- 
havior observed in impurity spin systems with nonuniform exchange couplings.- In fact, 
the magnetic nanostructures studied in this paper can be viewed as a super-spin glass 
with giant magnetic moments that interact via the dipole coupling. 

One of the experimental approaches to 2D nanostructures derives from investigations 
of three-dimensional (3D) ferrofluids consisting of non-overlapping magnetic particles in a 
liquid or frozen solvent.- The controlled preparation of 2D nanoparticle ensembles with dif- 
ferent degrees of structural and magnetic disorder has recently become feasible by reducing 
the thickness of 3D ferrofluids to within a few particle diameters.-i^i^ A second experimen- 
tal approach to 2D magnetic nanosystems stems from the study of thin magnetic films 
that are grown by depositing atoms on nonmagnetic substrates.-'^^'ii'i^ In the early stages 
of film growth, before a continuous layer is formed, it is possible to obtain a controlled 
distribution of magnetic islands that interact by long-range forces like dipole coupling or 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions in case of metallic substrates.— As 
a consequence of the preparation conditions, these magnetic nanostructures are usually 
characterized by the presence of a considerable degree of structural disorder. This mani- 
fests itself, for instance, in the particle-size dispersion or in the irregularity of the particle 
arrangement. 

The magnetic behavior of magnetic nanostructures is determined by the interplay 
between single-particle properties (e.g., magneto-crystalline or shape anisotropics) and 
interparticle effects (e.g., dipolar, RKKY or exchange interactions). For weakly cou- 
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pled systems the local contributions dominate, and the interactions can be handled as 
perturbations. In this case the magnetic relaxation is governed by single-particle fluc- 
tuations that can be described by a superparamagnetic Arrhenius-like model showing 
blocking effects at low temperatures.— However, in the most interesting case of strongly 
interacting magnetic nanostructures, for example, for densely packed nanoparticles, the 
single-particle approach is no longer appropriate.- Here the system exhibits a complicated, 
spin-glass-like behavior, characterized by a complex energy landscape, a large number of 
metastable minima, and a non-coUinear magnetization of the particle ensemble.— For 
these strongly interacting systems the magnetization dynamics is determined by a collec- 
tive response. In other words, changes in the magnetization direction of a given particle 
inevitably cause neighboring particles to vary their magnetic directions as well. During 
the transition from a local energy-minimum to another one, within a magnetic relaxation 
process, the magnetization directions of many particles change. Therefore, the collective 
behavior of the nanostructure has to be taken into account from the start. 

The disorder within the nanostructure and the anisotropic long-ranged nature of the 
dipole interaction represent a serious challenge for theoretical investigations so that simple 
analytical approaches are in general not applicable. Therefore, numerical simulations have 
been performed in order to achieve a detailed description of these systems.— Phenomeno- 
logical, macroscopic quantities like the field-cooled and zero-field-cooled magnetization 
curves, the linear and nonlinear susceptibilities, or the magnetic relaxation rate have 
been calculated.-*^*^'^ Moreover, the magnetic transitions between energy minima have 
been investigated by solving a continuum dynamical model including dissipation and ther- 
mal fluctuation terms (Langevin dynamics) as well as by using path-integral methods.— 
In order to take into account the diversity of nonequivalent trajectories, these calculations 
involve optimizations in a space of dimension d = nuk, where n is the number of particles, 
z/ the number of degrees of freedom, and k the number discretizations of the relaxation 
path. 

The main purpose of this paper is to investigate the collective microscopic processes 
occurring along the magnetic relaxations of disordered magnetic nanostructures. To this 
aim we focus on the elementary transitions along the minimum energy path (MEP) con- 
necting two neighboring local minima across a single saddle point. These fundamental 
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processes can be regarded as the microscopic analogy to the Barkhausen-jumps of the 
magnetization in hysteresis loops. The actual macroscopic relaxation is the result of a 
succession of such transitions between metastable states. Understanding them is therefore 
very important for the modelization of the long-time dynamics. One of the goals of this 
work is to analyze the properties of these magnetic rearrangements from a local perspec- 
tive by determining the variations A0j in the directions of the magnetic moments at each 
nanoparticle i. In particular a correlation range Ac is defined, that measures the spatial 
extent of the relevant collective processes as a function of experimentally important pa- 
rameters like particle density and degree of disorder. Furthermore, the distribution of the 
energy barriers is quantified in order to determine how they correlate with the changes in 
the magnetic configurations. 

The remainder of the paper is organized as follows. In Sec. [Ill the theoretical back- 
ground describing the interacting magnetic nanoparticle system is presented. The meth- 
ods for calculating saddle points, adjacent minima, relaxation paths and the resulting 
changes in the magnetic order along elementary transitions are outlined. In Sec. [Tlll rep- 
resentative results for the elementary transitions and the associated correlation ranges are 
reported and discussed. The distribution of the energy barriers, and their relation to the 
variation in the magnetic configurations are investigated. Finally, we conclude in Sec. IIVI 
by pointing out some of the potential implications of this study for magnetization reversal 
dynamics and high-density magnetic recording. 



II. MODEL AND CALCULATION METHOD 



of 
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The purpose of this section is to describe the main aspects of the present mode 
disordered magnetic nanostructures. Further details of the theory may be found in Ref. 
We consider a 2D system of non-overlapping, spherical magnetic particles contained 
in a square unit cell with periodic boundary conditions. Due to the strong interatomic 
ferromagnetic exchange interactions and the small size of the particles under consideration, 
it is clear that each particle i can be regarded as a single-domain or Stoner-Wohlfarth 
magnet carrying a giant spin nij.— The modulus of the local moments is given by rrii = 
ni fiat, where fiat is the atomic magnetic moment and rij the number of atoms in particle i. 
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Different kinds of structural arrangements are considered in the following: (i) a periodic 
square-lattice array, (ii) a disordered array where the particles are randomly dispersed 
around the sites of a periodic square lattice according to a Gaussian distribution with 
standard deviation ct/j, and (iii) a fully-random particle distribution of non-overlapping 
particles within the unit cell. 

In this paper we focus on the effects of the dipole coupling on the collective mag- 
netic behavior. This is a particularly interesting physical situation, since in the limit 
of weak magnetic anisotropy and external magnetic field the numbers of local minima 
and the complexity of the energy landscape are expected to be largest. For each particle 
arrangement the dipolar interaction between the particle moments nij is given by 



where Vij = — is the vector between the centers of particles i and j, rij = \rij\ 
the corresponding distance, and fio the vacuum permeability. The infinite range of the 
dipole interaction is taken into account by computing an Ewald-type summation over all 
periodically arranged unit cells of the extended planar system.— Since the particles do 
not have direct metal-metal contacts, no interparticle exchange interactions are included. 
Taking into account these short-range contributions, as well as local magnetic anisotropics 
or external magnetic fields, poses no major technical difficulties and could be investigated 
by similar simulations. For simplicity we assume that all particles have the same size 
Hi = 30000 atoms, which corresponds to a particle radius of about 15 lattice constants. 
Previous studies have actually shown that the dominant disorder effects are due to the 
irregularity of the particle arrangement and that the distribution of particle sizes yields 
a similar behavior.— The directions of the particle magnetizations irij are restricted to 
the plane of the nanostructure. Hence, the magnetic configurations are characterized for 
simplicity by the set of in-plane angles {(pi} in the unit cell (i = 1, . . . A^). 

The particle ensemble is characterized by the number of particles in the unit cell, 
by the standard deviation ct/j of the structural disorder of the particle array, and by 
the surface coverage C— The size of the unit cell is defined in terms of and the 
average interparticle distance Rq. For a given realization of the nanostructure we sample 
randomly the energy landscape and determine a large number of first-order saddle points 
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(SPs) or transition states. These are characterized by having one negative eigenvalue of 
the Hessian of the interaction energy E [see Eq. ([1])]. To locate the SPs we use a variant 
of the eigenvector-following method that is based on the eigenvalues and eigenvectors of 
the Hessian matrix as described in the Appendix. The two neighboring energy minima 
adjacent to each SP and the associated MEP are obtained by a simple steepest descent 
procedure starting at the SP along the negative-curvature direction. In this way the 
elementary relaxation path, the energies Es and Eui 2 for the saddle point S and adjacent 
minima Mi and M2, as well as the corresponding sets of magnetization angles {(j)i{S)} 
and {0i(Mi^2)}, are determined. In this context it is interesting to mention that loop-like 
paths also exist, where the SP connects a local minimum with itself via a closed loop 
(i.e.. Ml = M2). The presence of these features of the energy landscape is favored by the 
periodicity of the angular functions entering Eq. ([1]). Nevertheless, loop-like MEPs have 
been excluded from the present analysis, since they do not contribute to the magnetic 
relaxation. 

The range Ac of the magnetic correlations involved in an elementary transition is 
quantified in two steps as follows. First, the epicenter Rc of the angular variations 
A0j = |</)j(Mi) — </)j(M2)| between the two adjacent minima is determined from 

where rj is the position of the center of particle i. Note that Rc needs not to coincide 
with the position of any particle. Second, the dependence of A0j on the distance |rj — Rc| 
to the epicenter is fitted by an exponential of the form 

A^i = A^max exp(-2 |ri - Rd/Ac) , (3) 

from which the correlation range Ac is obtained. The angular variations are normalized 
with respect to the maximal change A0max occurring in that particular pair of states. The 
factor 2 in the exponential is introduced in order that Ac covers most of the significant 
angular variations. Hence, for |rj — Rd = Ac the change of angle A0j is reduced by a 
factor with respect to A^max- However, the use of an exponential fit should not be 
interpreted as a statement on the precise distance dependence of A</)j, which is presently 
not known in detail. 
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The Euclidean distance between the saddle point S and the adjacent minimum M is 
given by 



DsM provides a quantitative measure of the changes in the magnetic order during an 
elementary transition. The correlation ranges Ac and distances -Dsm are averaged over a 
large number of elementary transitions (typically a few thousands) and several geomet- 
rical arrangements of the particles in the nanostructure. In the following calculations 
parameters appropriate for Fe particles are used, i.e., fiat = 2.2 /i^ and nearest neighbor 
distance Cq = 2.5 A. 

III. RESULTS AND DISCUSSION 

The most stable magnetic arrangement of a periodic square lattice of dipole coupled 
magnetic nanoparticles having all the same size is know to correspond to the so-called 
microvortex state.— This magnetic order has a vanishing net magnetization and is charac- 
terized by a continuous degeneracy with respect to the microvortex angle. Any, however 
small, deviation from the square-lattice symmetry lifts this degeneracy. For weak disorder 
the low-energy states preserve a close resemblance with a microvortex state but as the dis- 
order increases the apparent extended periodicity of the magnetic structure is completely 
lost. In the case of strong disorder the magnetic state is characterized by head-to-tail 
arrangements of the particle moments and by local flux closures of the magnetization 
profile. In this Section results are presented and discussed for the correlation range Ac, 
the energy barriers AE^m = — Em, and the distances -Dsm, including averages over 
many elementary transitions as well as different values of the parameters A^, and C 
characterizing the magnetic nanostructure. 

Figure [T] illustrates typical examples of elementary magnetic transitions. The arrows 
refer to the magnetic directions of two adjacent minima, and the different shadings in- 
dicate the importance of the associated angular changes A0j. Two different structural 
arrangements are considered: a) an almost periodic square-lattice array with = 5%, 
and b) a fully-random particle setup. The resulting correlation range Ac is indicated by a 
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FIG. 1: (Color online) Illustration of elementary transitions in dipole-coupled nanostructures 
having a) a weak disorder with an = 5%, and b) a random particle distribution. In both cases 
the unit cell (dotted frame) contains N = 100 particles and the coverage is C = 0.35. The arrows 
show the directions of the magnetic moments corresponding to two adjacent local minima that 
are separated by a single energy barrier. The shading indicates the importance of the changes 
in the direction of the local moments. The dashed circle refers to the correlation range Ac where 
the angular variations are most significant. 
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dashed circle. One clearly observes that Ac is significantly larger for the almost periodic 
arrangement than for the random ensemble. In fact, in the former case Ac is often compa- 
rable with the size of the unit cell. Hence, to obtain reliable quantitative results one would 
have to consider a much larger cell so that spurious effects from the periodic boundary 
conditions are avoided. In contrast, for strong disorder the region with significant orien- 
tational changes of rrij remains quite small. In this case the results for Ac are not affected 
by considering even larger unit cells. The magnetic transition presented in Figure Wljo) 
refers to a relatively large Ac, usually the regions of significant magnetic variations for this 
kind of arrangement are somewhat smaller. Note that for small disorder. Figure [T](a), the 
magnetic variations are not isotropic but still refiect the square symmetry. 

In order to derive a quantitative estimate of Ac we approximate the dependence of A</)j 
on |rj — Rc| by an exponential law [see Eq. ([3])]. Figure [2] shows the correlation range 
Ac averaged over a large number of elementary transitions as a function of the structural 
disorder an. A single realization of the nanostructure is considered either for different 
coverages C keeping the number of particles in the unit cell constant, or vice versa. 



Results are given for Ac and for the standard deviation AAc = y ((Ac — Ac)^) in units 
of the average interparticle distance Rq. Notice that Rq increases as the coverage C is 
reduced. One observes that both Ac and AAc decrease rapidly with increasing ct/j. Already 
for cT/j > 10-20% their values are practically indistinguishable from those of a random 
particle ensemble. Usually Xc/Ro does not depend strongly on the coverage C, except 
for the largest values (e.g., C = 0.55) where the non-overlap constraint and the assumed 
small size dispersion force the particles to adopt an almost periodic arrangement. All the 
calculations for different numbers A^ of particles in the unit cell yield very similar results 
within the statistical uncertainty. In particular the elementary transitions illustrated in 
Fig. [1] correspond to relatively large correlation ranges. Moreover, notice that for small 
structural disorder ctr, < 5% and for small A^ our results for Ac might not be very accurate 
quantitatively, since in these cases the correlation range is comparable to the size of the 
unit cell. 

The regions where the elementary transitions take place become more localized with 
increasing disorder of the nanostructure. Beyond Ac, in the remaining part of the unit 
cell, no significant changes in the directions of the particle magnetizations are found. The 




9 





20 30 

C^R [%] 



FIG. 2: (Color online) Average correlation range Ac and standard deviation AAc as functions of 
the structural disorder aR. The results are obtained by sampling a large number of elementary 
transitions. In the top figure different coverages C are considered, while in the bottom figure 
the number of particles N in the unit cell is varied. 

magnetic transitions in strongly interacting particle ensembles can be studied by using 
periodic cells provided that they are significantly larger than Ac. For strong disorder al- 
ready small unit cells suffice, whereas for almost periodic ensembles much larger cells are 
required. This renders the simulations far more demanding since an important number 
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FIG. 3: (Color online) Distribution of the energy barriers AE^m corresponding to the elemen- 
tary transitions in a random ensemble having a coverage C = 0.35. The curve shows a fitted 
exponential law. 

of particles takes part in an elementary transition. In the limit of a periodic arrangement 
the whole system is involved in an elementary relaxation (microvortex state). Neverthe- 
less, for moderately disturbed ensembles where the underlying square arrangement is still 
clearly recognizable, the correlation range already assumes relatively small values [see, 
for example. Fig. [I](a)]. It would be therefore interesting to investigate the dependence 
of Ac on an in more detail, particularly as the periodic lattice is approached {api 0). 
Finally, it should be mentioned that including single-particle magnetic anisotropies in the 
calculations would tend to further reduce Ac, especially for small coverages C and large 
interparticle distances Rq. This corresponds to weakly interacting particle ensembles, 
where the elementary transitions are essentially given by the magnetic reversal of a single 
particle. 

The energy barriers AEsu = Es — Eu between a saddle point and the adjacent minima 
are one of the central properties governing the relaxation dynamics. In Fig. [3] the distri- 
bution of AEsM is shown for a random ensemble having a coverage C = 0.35. The results 
were derived taking into account 826 different elementary transitions of the same nanos- 
tructure. One observes that the probability density of finding an energy barrier AEsm 
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FIG. 4: (Color online) Scatter plot of the energy barrier AE^sM as a function of the distance 
DsM between a saddle point and the adjacent local minima. As in Fig. [3] the particle setup is 
random and the coverage is C = 0.35. 

follows approximately an exponential behavior, although the precise energy dependence 
of the distribution is not known. In any case, it is worth noting that the distribution 
of barrier energies contrasts with those of the absolute energies of the saddle points and 
local minima, which were shown to be Gaussian-like for random particle ensembles.— 

Figure m shows the correlation between AE^u and the distance -Dsm between the saddle 
points and the adjacent minima in a random nanoparticle setup [see Eq.(jl])]. As a general 
trend one observes that if D^m is small also the corresponding AE^m is small. However, for 
larger -Dsm the energy barriers can assume both large and small values. Hence, important 
changes of the magnetic arrangement are also possible by involving only a single and 
relatively small energy barrier AEsu- In this connection it should be noted that the 
distribution of the distances -Dsm between saddles and minima has been found to follow 
a Poisson-like distribution.— 

In Fig. \5\ the average distance -Dsm and the standard deviation A-Dsm = 



Y ((-Dsm — -Dsm)^) are given as functions of the structural disorder aR for different cover- 
ages C. -Dsm provides a measure of the change in the magnetic state associated with an 
elementary transition. In agreement with the results for Ac shown in Fig. [2], -Dsm decreases 
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FIG. 5: (Color online) Average distance Dsm between saddle points and adjacent minima, 



as well as the standard deviation ADsm = y ((-^sm — -^sm)^) as functions of the structural 
disorder cr/j. The results are derived from a large number of elementary transitions for each 
different coverage C. 

with increasing a^i, converging to the random-setup hmit aheady for ct/j > 20 %. No sig- 
nificant dependence on the coverage C is found. The distance between the neighboring 
minima of elementary transitions (-Dmm) follows closely the behavior of -Dsm? of course 
with somewhat larger values. As in the case of Fig. [2] the results for weak structural 
disorder < 5% are less reliable due to the finite size of the periodic cell. 

IV. CONCLUSION 

The elementary relaxation processes occurring in disordered 2D ensembles of dipole- 
coupled magnetic nanoparticles have been investigated by calculating the transitions be- 
tween neighboring metastable states across first-order saddle points. For this purpose an 
eigenvector-following method was implemented. The extension of the correlated region 
where the magnetic directions of the particles vary significantly has been calculated. In 
particular we have shown that the magnetic correlation range Ac becomes smaller as the 
degree of disorder in the nanostructure increases. The tendency to localization of elemen- 
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tary transitions in disordered magnetic systems renders the simulation of these structures 
computationally less demanding. Moreover, the distribution of the energy barriers AEsu 
has been found to follow approximately an exponential behavior. An analysis of the rela- 
tion between Ai^sM and the change -Dsm of the magnetic state shows that small energy 
barriers can lead to both small and large magnetic relaxations, whereas large barriers 
always imply an important rearrangement of the magnetic state. 

A number of implications and open questions related to our results may be pointed out. 
As already noted, the time evolution of the magnetization of the nanostructure can be 
described by a succession of the elementary transitions investigated in this work. Partic- 
ularly important for the description of such processes is the knowledge of the connectivity 
of the network of metastable states in the complex spin-glass-like energy landscape. For 
each given system one would have to determine which minima are connected by elemen- 
tary transitions, as well as the corresponding transition rates. The latter are mainly 
defined by the energy barrier, the eigenfrequencies at the saddle point, and the length 
of the MEP. They can be determined in the framework of transition-state theory in- 
cluding finite-temperature effects. Besides the elementary transitions involving collective 
rotations of the particles moments within a limited correlation range, the magnetic re- 
laxation of the nanostructure as a whole involves some complex domain-wall-like motion 
that results from a sequence of elementary processes. In the case of time-dependent fields, 
averages over trajectories should be performed in order to simulate hysteresis loops from 
a microscopic point of view. 

For more realistic simulations and for a detailed comparison with experiment it would 
be interesting to take into account other potentially important contributions, for ex- 
ample, the magnetic anisotropy energy of individual particles and the interaction with 
external magnetic fields. Qualitatively, one expects that increasing the single-particle 
anisotropies or the applied field should tend to reduce the correlation range of the elemen- 
tary dynamical process, since the magnetic relaxation would be increasingly dominated by 
single-particle effects. Another aspect deserving further study concerns the distribution 
of particle sizes, which is certainly unavoidable in real nanostructures. In fact previous 
works2i suggest that they should have similar, probably less severe consequences on the 
energy landscape as the positional disorder considered here. More interesting interaction 
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effects on the dynamics are expected to result from other types of interparticle couphngs, 
for instance, indirect RKKY interactions between NP on metaUic substrates, or direct ex- 
change interactions at the boundaries of bare NPs in contact. These energy contributions 
can be treated straightforwardly using the methods developed in this work. 

The trend to increasing localization of the magnetic excitations with increasing dis- 
order in the particle arrangement should be of interest for tailoring materials for infor- 
mation storage. In fact, a stronger disorder decreases the minimum distance required for 
an independent switching of neighboring 'bits'. Hence, controlling the main parameters 
characterizing the nanostructure allows to optimize their magnetic response as potential 
high-density magnetic-storage media or spintronic devices, for example, by tuning the par- 
ticle arrangements and the coverages. This also reveals the limits for single-particle read 
and write processes in interacting magnetic nanostructures. Research in these directions 
is currently in progress. 



APPENDIX A: ITERATIVE SEARCH FOR SADDLE POINTS 



The method to determine the first-order saddle points is based on the eigenvalues and 
eigenvectors of the Hessian matrix of the energy E a.s a. function of the magnetization 
angles [see Eq. ([T])].— The following iteration algorithm is applied as sketched in 
Fig. EJ (i) Consider an arbitrary initial magnetic state P^'^^ in the A^- dimensional space of 
the magnetization angles (ii) Determine at P^'""-' the eigenvector Vq''"'' of the Hessian 

corresponding to the lowest eigenvalue Xl^\ (iii) Maximize the energy of the magnetic 
state P along Vq'^'* until a state Pq is found where the energy shows a maximum along 
Vq'^'*, or until a maximal prescribed distance from P^''"^ is reached, (iv) Starting now from 
Pq minimize the energy in the (A^ — l)-dimensional space perpendicular to Vq'^'* until the 
nearest minimum in this restricted subspace is reached, (v) Set this minimum as the initial 
state p('^+^) of the next iteration step, and repeat steps (ii)-(iv) until a saddle point S is 
reached with the desired accuracy. This procedure is similar to the eigenvector-following 



method described in Ref. 
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The convergence criterion at a first-order saddle point S requires that the norm of the 
gradient be sufficiently small, that the lowest eigenvalue Aq at S be negative, and that 
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FIG. 6: (Color online) Sketch of the eigenvector-following method to determine first-order 
saddle points. 

all other eigenvalues be positive. The two corresponding minima Mi and M2 adjacent to 
a saddle point are then determined by descending from the saddle point along the two 
opposite directions defined by the eigenvector vq corresponding to Aq. The resulting MEP 
is obtained by following the gradient by steepest descent. 
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